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Abstract 

Cylindrically symmetric inhomogeneous cosmological model for bulk 
viscous fluid distribution with electromagnetic field is obtained. The 
source of the magnetic field is due to an electric current produced along 
the z-axis. F12 is the non-vanishing component of electromagnetic field 
tensor. To get the deterministic solution, it has been assumed that the 
expansion 9 in the model is proportional to the shear a. The values of 
cosmological constant for these models are found to be small and positive 
at late time which are consistent with the results from recent supernovae 
la observations. Physical and geometric aspects of the models are also 
discussed in presence and absence of magnetic field. 

Keywords : cosmology, variable cosmological term, electromagnetic field, inho- 
mogeneous universe 
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1 Introduction 

Inhomogeneous cosmological models play an important role in understanding 
some essential features of the universe such as the formation of galaxies during 
the early stages of evolution and process of homogenization. The early attempts 
at the construction of such models have been done by Tolman [1] and Bondi [2] 
who considered spherically symmetric models. Inhomogeneous plane-symmetric 
models were considered by Taub [3J H] and later by Tomimura [5] , Szekeres [B] , 
Collins and Szafron |7| , Szafron and Collins [8] . Recently, Senovilla [9] obtained 
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a new class of exact solutions of Einstein's equations without big bang singular- 
ity, representing a cylindrically symmetric, inhomogeneous cosmological model 
filled with perfect fluid which is smooth and regular everywhere satisfying en- 
ergy and causality conditions. Later, Ruiz and Senovilla [10] have examined 
a fairly large class of singularity free models through a comprehensive study 
of general cylindrically symmetric metric with separable function of r and t as 
metric coefficients. Dadhich et al. [11] have established a link between the FRW 
model and the singularity free family by deducing the latter through a natural 
and simple in-homogcnization and anisotropization of the former. Also, Patel 
ct al. [fj] presented a general class of inhomogeneous cosmological models filled 
with non-thermalized perfect fluid by assuming that the background space-time 
admits two space-like commuting Killing vectors and has separable metric co- 
efficients. Singh, Mehta and Gupta [13] obtained inhomogeneous cosmological 
models of perfect fluid distribution with electro-magnetic field. Recently, Prad- 
han et al. [14] have investigated plane-symmetric inhomogeneous cosmological 
models in various contexts. Cylindrically symmetric space-time play an im- 
portant role in the study of the universe on a scale in which anisotropy and 
inhomogeneity are not ignored. Roy and Singh [15] , Bali and Tyagi [16[ 117]. 
Chakrabarty et al. [T8] and Pradhan et al. [TI5] have investigated cylindrically 
symmetric inhomogeneous cosmological models in presence of electromagnetic 
field. 

The occurrence of magnetic field on galactic scale is well-established fact 
today, and their importance for a variety of astrophysical phenomena is gener- 
ally acknowledged as pointed out by Zeldovich et al. [20]. Also Harrison [21] 
has suggested that magnetic field could have a cosmological origin. As a nat- 
ural consequences, we should include magnetic fields in the energy-momentum 
tensor of the early universe. The choice of anisotropic cosmological models in 
Einstein system of field equations leads to the cosmological models more gen- 
eral than Robertson- Walker model [22]- The presence of primordial magnetic 
field in the early stages of the evolution of the universe has been discussed by 
several authors [53] — [35]. Strong magnetic field can be created due to adia- 
batic compression in clusters of galaxies. Large-scale magnetic field gives rise to 
anisotropies in the universe. The anisotropic pressure created by the magnetic 
fields dominates the evolution of the shear anisotropy and it decays slower than 
the case when the pressure was isotropic [331 131] ■ Such fields can be gener- 
ated at the end of an inflationary epoch [35] — [39]. Anisotropic magnetic field 
models have significant contribution in the evolution of galaxies and stellar ob- 
jects. Bali and Ali [40j obtained a magnetized cylindrically symmetric universe 
with an electrically neutral perfect fluid as the source of matter. Chakrabarty 
et al. ,Tj5j and Pradhan et al. [H] have investigated magnetized viscous fluid 
cosmological models in various contexts. 

There are significant observational evidence for the detection of Einstein's 
cosmological constant, A or a component of material content of the universe that 
varies slowly with time and space to act like A. Some of the recent discussions 
on the cosmological constant "problem" and on cosmology with a time-varying 
cosmological constant by Ratra and Peebles [42] , and Sahni and Starobinsky 
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[13] point out that in the absence of any interaction with matter or radiation, 
the cosmological constant remains a "constant". However, in the presence of 
interactions with matter or radiation, a solution of Einstein equations and the 
assumed equation of covariant conservation of stress-energy with a time- varying 
A can be found. This entails that energy has to be conserved by a decrease in the 
energy density of the vacuum component followed by a corresponding increase 
in the energy density of matter or radiation (see also Carroll, Press and Turner 
[H] , Peebles [JS] , Padmanabhan [46 ) . There is a plethora of astrophysical ev- 
idence today, from supernovae measurements (Perlmutter et al. [17], Riess et 
al. [IS], Garnavich et al. [49] . Schmidt et al. [50], Blakeslee et al. [HI], Astier 
et al. [52]), the spectrum of fluctuations in the Cosmic Microwave Background 
(CMB) [53], baryon oscillations [54] and other astrophysical data, indicating 
that the expansion of the universe is currently accelerating. The energy budget 
of the universe seems to be dominated at the present epoch by a mysterious dark 
energy component, but the precise nature of this energy is still unknown. Many 
theoretical models provide possible explanations for the dark energy, ranging 
from a cosmological term |55j to super-horizon perturbations |56] and time- 
varying quintessence scenarios [57] . These recent observations strongly favour 
a significant and a positive value of A with magnitude A(Gh/c 3 ) rs 1CP 123 . In 
Ref. [15J, Riess et al. have recently presented an analysis of 156 SNe includ- 
ing a few at z > 1.3 from the Hubble Space Telescope (HST) "GOOD ACS" 
Treasury survey. They conclude to the evidence for present acceleration qo < 
(q « —0.7). Observations (Knop et al. [58] ; Riess et al., [48]) of Type la Su- 
pernovae (SNe) allow us to probe the expansion history of the universe leading 
to the conclusion that the expansion of the universe is accelerating. 

Most studies in cosmology involve a perfect fluid. Large entropy per baryon 
and the remarkable degree of isotropy of the cosmic microwave background radi- 
ation, suggest that we should analyze dissipative effects in cosmology. Further, 
there are several processes which are expected to give rise to viscous effect. 
These are the decoupling of neutrinos during the radiation era and the recom- 
bination era [50J , decay of massive super string modes into massless modes [BOJ , 
gravitational string production [ST] [02] and particle creation effect in grand uni- 
fication era [63] . It is known that the introduction of bulk viscosity can avoid 
the big bang singularity. Thus, we should consider the presence of a material 
distribution other than a perfect fluid to have realistic cosmological models (see 
Gr0n [64] for a review on cosmological models with bulk viscosity). A uniform 
cosmological model filled with fluid which possesses pressure and second (bulk) 
viscosity was developed by Murphy [65] . The solutions that he found exhibit 
an interesting feature that the big bang type singularity appears in the infinite 
past. 

Recently, Pradhan et al. |66] have obtained inhomogeneous perfect fluid 
universe with electromagnetic field. Motivated by the situation discussed above, 
in this paper, we have obtained a new cylindrically symmetric inhomogeneous 
cosmological model for bulk viscous fluid distribution in presence and absence of 
electromagnetic field. The coefficient of bulk viscosity is assumed to be a power 
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function of mass density. This work generalize the previous work of Pradhan 
et al. |66j . Some physical and geometric behaviour of the models in presence 
and absence of magnetic field are also discussed. This paper is organized as 
follows. The metric and the field equations are laid down in Section 2. In 
Section 3, we deal with the solutions of the field equations in presence of bulk 
viscous fluid with electromagnetic field and variable cosmological term. We have 
also described the physical and geometric aspects of the models. In Section 4, 
we obtain the solutions of the field equations in absence of the magnetic field. 
Finally in Section 5 concluding remarks are given. 

2 The Metric and Field Equations 

We consider the metric in the form 

ds 2 =A 2 {dx 2 -dt 2 ) + B 2 dy 2 + C 2 dz 2 , (1) 

where A is the function of t alone and B and C are functions of x and t. The 
energy momentum tensor is taken as has the form 

if = (P + P)uiU 3 + pg{ + Ef, (2) 

where Ef is the electromagnetic field given by Lichnerowicz [67j 

h x ti L iU j + hA - hih? , (3) 



and 



P = P-£,<- (4) 



Here p, p, p and £ are the energy density, isotropic pressure, effective pressure 
and bulk viscous coefficient respectively and u l is the fluid four-velocity vector 
satisfying the condition 

giju'u 1 = -1. (5) 
p is the magnetic permeability and hi, the magnetic flux vector defined by 

hi = ^*F jiU j , (6) 
P 

where the dual electromagnetic field tensor *Fy is defined by Synge [68] 



= —^-£i 3 kiF kl . (7) 

Here F^ is the electromagnetic field tensor and tijki is the Levi-Civita tensor 
density. 

The co-ordinates are considered to be comoving so that u 1 = u 2 = u 3 = and 
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u 4 = -4. If we consider that the current flows along the z-axis, then F12 is the 
only non- vanishing component of Fij . The Maxwell's equations 

Fij-k + Fjk-i + Fki-j = 0, (8) 
= 4tt ,F (9) 

■3 

require that F12 is the function of x-alone. We assume that the magnetic per- 
meability is the functions of x and t both. Here the semicolon represents a 
covariant differentiation. 



The Einstein's field equations (in gravitational units c — 1, G = 1) read 



1 



Rl - -Rffl + Ag? = -87tf* 
for the line-element (JXJ) lead to the following system of equations: 
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where the sub indices 1 and 4 in A, B, C and elsewhere denote ordinary differ- 
entiation with respect to x and t respectively. 



3 Solution of the Field Equations 

Equations (| 1 1 1) - fj 1 5 j) are five independent equations in six unknowns A, B, C, 
p, p and Fi2- For the complete determinacy of the system, we need one extra 
condition. The research on exact solutions is based on some physically reason- 
able restrictions used to simplify the Einstein equations. 
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To get determinate solution we assume that the expansion 6 in the model is 
proportional to the shear a. This condition leads to 

A=(^)\ (16) 



where n is a constant. From Eqs. pip - (fT5|) . we have 

d n 

K (constant) (17) 
(18) 

We also assume that 

B = f(x)g(t) 

C = f(x)k(t). (19) 
Using Eqs. JIB]) and ([17]) in ([15]) and dTTJ) lead to 
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ffu - fx = ^/ 2 - (22) 

Equation (|20ll leads to 

k = cg a , (23) 

where a = 2n+i ano - c * s ^ ne cons t an t of integration. From Eqs. (f2"Tj) and (f2"3"]) , 

we have 

^+/4 = Af, (24) 

where 

na(a — 1) + a K 
1 7 T~\ i -I ' 



n(a- 1) + 1 ' n(l - a) - 1 

Equation (|22|) leads to 

T 



/ = cxp\^K(x + x ) 2 ), (25) 
where Xo is an integrating constant. Equation (f2"4")l leads to 

.9= ( Cl e b * + C2 e- W )^, (26) 

where 6 = \fTj3 + 1)N and ci, c 2 are integrating constants. Hence from (j2"3"j) 
and (|26|) . we have 

fc = c(cie bt +c 2 e- bt )"* T ) . (27) 
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Therefore we obtain 
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where a — — , C3 being a constant of integration. 

After using suitable transformation of the co-ordinates, the model (JXJ) reduces 
to the form 

ds 2 = a 2 ( Cl e bT + c 2 e- 6T )W(rfI 2 - dT 2 ) + e KX \c x e bT + c 2 e- bT )W^dY 2 



+ e K * ( Cl e bT + c 2 e- bl ) — dZ 



(31) 



where x + xq = X, t = T, y = Y, cz = Z. 

The expressions for effective pressure p and density p for the model (I3ip are 
given by 
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For the specification of £, we assume that the fluid obeys an equation of state 
of the form 

P = IP, (34) 

where 7(0 < 7 < 1) is a constant. Thus, given £(t) we can solve for the 
cosmological parameters. In most of the investigation involving bulk viscosity 
it is assumed to be a simple power function of the energy density [69) — [73] 



(35) 



where £0 and m are constants. For small density, m may even be equal to unity 
as used in Murphy's work [65) for simplicity. If m = 1, (]35|) may correspond 
to a radiative fluid [75] . Near a big bang, < m < h is a more appropriate 
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assumption [TJ] to obtain realistic models. 
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where 9 is the scalar of expansion calculated for the flow vector u l and is given 
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For simplicity and realistic models of physical importance, we consider the fol 
lowing two cases (m = 0, 1). On using ([33|) in (f3"2"|) , we obtain 



3.1 Model I: Solution when m = 
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Figure 2: The plot of cosmological term A(T) Vs. T 



3.2 Model II: Solution when m = 1 

When m = 1, Eq. (|55j) reduces to £ = £ p. With the use of Eqs. 
103), Eq. ([21]) reduces to 
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From Eqs. (|39f and (|4"Tj) , we note that is a decreasing function of time 
and p > for all times. This behaviour is clearly depicted in Figures 1 and 3 
as a representative case with appropriate choice of constants of integration and 
other physical parameters using reasonably well known situations. Figures 1 
and 3 show this physical behaviours of energy density as a decreasing functions 
of coordinate time T of Model I and II. Here the coordinate time T is related to 
normal time as shown in Ref. [19] ■ This also follows for rest part of the paper. 



In spite of homogeneity at large scale our universe is inhomogeneous at small 
scales, so physical quantities being position dependent are more natural in our 
observable universe if we do not go to super high scale. This result shows this 
kind of physical importance. In recent time the A-term has interested theoreti- 
cians and observers for various reasons. The nontrivial role of the vacuum in the 
early universe generate a A-term that leads to inflationary phase. Observation- 
ally, this term provides an additional parameter to accommodate conflicting 
data on the values of the Hubble constant, the deceleration parameter, the 
density parameter and the age of the universe (for example, see the references 
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Figure 3: The plot of energy density p(T) Vs. T 




Figure 4: The plot of cosmological term A(T) Vs. T 
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[731(73]). Assuming that A owes its origin to vacuum interactions, as suggested 
in particular by Sakharov [77], it follows that it would in general be a function 
of space and time coordinates, rather than a strict constant. In a homogeneous 
universe A will be at most time dependent [78j . In our case this approach can 
generate A that varies both with space and time. In considering the nature of 
local massive objects, however, the space dependence of A cannot be ignored. 
For details discussion, the readers are advised to see the references (Narlikar, 
Pecker and Vigier |79J, Ray and Ray 80J, Tiwari, Ray and Bhadra [8"T]). 

The behaviour of the universe in this model will be determined by the cos- 
mological term A ; this term has the same effect as a uniform mass density 
p e ff — —A/4ttG, which is constant in space and time. A positive value of A 
corresponds to a negative effective mass density (repulsion). Hence, we expect 
that in the universe with a positive value of A, the expansion will tend to ac- 
celerate; whereas in the universe with negative value of A, the expansion will 
slow down, stop and reverse. From Eqs. (I40[) and (U2J), we see that the cos- 
mological term A is a decreasing function of time and it approaches a small 
positive value at late time. From Figures 2 and 4, we note this behaviour of 
cosmological term A in both models I and II. Recent cosmological observations 
suggest the existence of a positive cosmological constant A with the magnitude 
A(Gh/c 3 ) ~ 10 -123 . These observations on magnitude and red-shift of type la 
supernova suggest that our universe may be an accelerating one with induced 
cosmological density through the cosmological A-term. Thus, our model is con- 
sistent with the results of recent observations. 



Some Physical and Geometric Features : 

The non-vanishing component F12 of electromagnetic field tensor is obtained as 



Jx 6 2 (1 - a) 



,ffX' 



8tt 09 + 1) 

(fi-a) ( Cl e bT 



( Cl e w + C2 e- w )WiTx 



1 



c 2 e-» T ) 2 



((3+1) ( Cl e bT + c 2 e- bT ) 2 



(43) 



The expressions for the shear scalar a 2 , acceleration vector Uj and proper volume 
V 3 for model ([31]) are given by 
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The rotation u> is identically zero. From set of equations (|3"T|) , (|4"3"| - (|4T|) . the 
model brings out the following features: 

The model starts expanding at T > and goes on expanding indefinitely when 
7 \p+i) ^ 0' m °del pip represents an expanding, shearing and non-rotating 
universe in which the flow vector is geodetic. Since ? = constant, the model 
does not approach isotropy. As T increases the proper volume also increases. 
The model is non-accelerating. The physical quantities p and p decrease as F\2 
increases. However, if Tg^jy > 0, the process of contraction starts at T > 
and at T — oo the expansion stops. The electromagnetic field tensor does not 
vanish when 6^0, and 1. 



4 Solution in Absence of Magnetic Field 

In absence of magnetic field the Einstein field equations for metric (JlJ read as 
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which on integration gives 



g = (c 4 t + c 5 ) <=+d , 
where C4 and C5 are constants of integration. Hence from 

k = c(at + C 5 )T^TT) _ 
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In this case (|T71) leads to 
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Therefore, we have 
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where a is already defined in previous section. 

After using suitable transformation of the co-ordinates, the model ([T]) reduces 
to the form 
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The expressions for effective pressure p and density p for the model (|6"2"|) are 
given by 
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n(l — a 2 ) + a 1 
(a + 1) 2 T 2 



i\T(2 + 3KX 2 ) 



On using (|35|) in (|63|) , we obtain 
87r(p-£ o p m 0) = J 



2ti(1-q) 



a 2 (c 4 T)" 



n(l-a 2 ) + a 1 2 2 ' 
(a + l)2 J^+ KX 



A, (63) 
+ A. (64) 

A, (65) 



where 9, in this case, is calculated for the flow vector u 1 and is given by 



where 



„(l- Q ) + (l + c) ' 

T 



n(l - a) + (1 + a) 



(66) 
(67) 



a(l + a)c 4 ( 
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4.1 Model I: Solution when m = 



When m = 0, Eq. (j3"5j) reduces to £ = £ - With the use of Eqs. 
(|6"6"j) , Eq. IBS]) reduces to 



and 



4tt(1+7)p 



1 



2n(l-a) 

a 2 (c 4 T) (!+«.)' 



n(l — a 2 ) + a 1 



(a + 1) 2 T 2 
Eliminating between ([Ml) and (f6"8")) . we get 



- iT(l + ifX 2 ) 



Air to K 3 



■»(l- a ) + (l + oQ 

j 1 fiTHj 



(68) 



Figure 5: The plot of energy density p(T) Vs. T 



(1+ 7 )A 



1 



2re(l-a) 



a 2 (c 4 T)- 



{n(l-a 2 ) + a} (l- 7 ) 



y2 



8ir£ K 3 

n(l-c) + (l + c) 



if{KI 2 (l + 37) + 27} 



(69) 



From Eq. |68|) . we see that p(f) is a decreasing function of time and p > for 
all times. Figure 5 shows this behaviour of energy density in Model I. 
From Eq. (f69|) . we observe that the cosmological term A is a decreasing function 
of time and it approaches a small positive value at late time. From Figure 6, 
we note this behaviour of cosmological term A in both model I. 
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Figure 6: The plot of cosmological term A(T) Vs. T 



4.2 Model II: Solution when m — 1 

When m = 1, Eq. (|35[) reduces to £ = £oP- With the use of Eqs. 
(f66|) , Eq. (|65j) reduces to 

i 

47T/3 = 



K 3 



„(i- Q) + (i +t> ) 



ri(l — a 2 ) + a 1 



(a + 1) 2 T 2 
Eliminating p(i) between ([Ml) and (|70)l . we get 

2 



- iT(l + i^X 2 ) 



A = 



a 2 (c 4 T)^w- l+ 7 -Co 



A':, 



i(l-c.) + (l + a ) 
T (! + °) 



?i(l - a 2 ) + a 1 



(a + 1) 2 T 2 

n(l — a 2 ) + a 1 



2~(l-n) 

a 2 (c 4 T) <i+°) 



-if(l + KI 2 ) 

K(2 + 3JOT 2 ) 



(a + 1) 2 T 2 



and 



(70) 



(71) 



From Eq. (|70|) and also from Figure 7(a), it seems that in very early stage of 
universe the energy density could be negative and may be link with some of the 



1G 



4 6 8 

T -» 



(a) 



(b) 



4 6 8 10 

T -» 



Figure 7: The plot of (a) p as a function of T, (b) A as a function of T 



early universe physics which is the domain of quantum cosmology or early quan- 
tum mechanics. So we do not infer any things about energy density. But note 
that once energy density is negative in the initial stage even after oscillations it 
returns back to small negative value at the later stage of the evolution. Hence 
the details of macro physics in early universe decides its fate. So we cannot 
infer from this model. Similar behaviours is also reflected about the cosmolog- 
ical constant A. We do not know for wide range of parameters this behaviour 
holds. It may be feasible that Figures 7(a,b) may coincide with Figures 3 and 
4. But it does not seem to be vary generic behaviour of generalized model of 
section 3.2. So we do not go to any conclusion about its applicability in general. 
But it seems that magnetic field and bulk viscosity prevent this kind of unusual 
behaviours as it is clear from Figures 1 to 4. 

Some Physical and Geometric Features: 
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The expressions for the shear scalar a 2 , acceleration vector iii and proper 
volume V 3 for model ([6"2")1 are given by 



2 _ {n(l - a) + (1 + a)} 2 - 3n(l - a 2 ) - 3a 

2ra(l-o) an(l-a)+2(l+a) 

3(l + a) 2 a 2 c 4 (1+Q> T 



6» 2 3{n(l-a) + (l + a)} 2 

The rotation u> is identically zero. 



(72) 



u i = (0,0,0,0), (73) 

^ 3 = = a 2 e KX (c 4 T) ~ . (74) 

From dZH) and ([66]). we obtain 

cr 2 _ {n(l - a) + (1 + a)} 2 - 3n(l - a 2 ) - 3a 



(75) 



The model in absence of magnetic field starts expanding with a big bang at 
T = and it stops expanding at T = oo. In absence of magnetic field, the model 
in general represents an expanding, shearing and non-rotating in which the flow 
vector is geodetic. Since 5 = constant, the model does not approach isotropy. 
As T increases the proper volume also increases. The model is non-accelerating. 



5 Concluding Remarks 

We have obtained a new cylindrically symmetric inhomogeneous cosmological 
model of electro-magnetic bulk viscous fluid as the source of matter. Generally 
the model represents expanding, shearing and non- rotating universe in which 
the flow vector is geodetic. It is worth mention here that in presence of mag- 
netic field the model (|3"Tj) is expanding whereas in absence of magnetic field the 
model ([6"2"|) starts with a big bang singularity. In both models § = constant 
and hence they do not approach isotropy. The models are non-accelerating. In 
these solutions all physical quantities depend on at most one space co-ordinate 
and time. 

It is important to note here that both the models (|31|) and ([62]) in presence and 
absence of magnetic field reduce to homogeneous universe when K = 0. This 
shows that for K — 0, inhomogeneity dies out. 

The effect of bulk viscosity is to produce a change in perfect fluid and hence 
exhibit essential influence on the character of the solution. The effect is clearly 
visible on the p effective (see details in previous sections). In Section 3, we 
have shown regular well behaviour of energy density, cosmological constant (A) 
and the expansion of the universe with parameter T. The section 4 is a toy 
investigation to see that the effect of bulk viscosity plays dynamic role in the 
evolution equations. We also observe that Murphy's conclusion [65] about the 
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absence of a big bang type singularity in the infinite past in models with bulk 
viscous fluid, in general, is not true. The results obtained by Myung and Cho 
[BU] also show that, it is, in general, not valid, since for some cases big bang 
singularity occurs in finite past. 

In presence and absence of magnetic field, the cosmological terms in models 
are decreasing function of time and approach a small value at late time (with 
exception to = 1 in the absence of magnetic field) . The values of cosmological 
"constant" for the models are found to be small and positive, as obtained in 
recent results from the supernovae observations recently obtained by the High- 
Z supernovae la Team and Supernovae Cosmological Project. Our solutions 
generalize the solutions recently obtained by Pradhan et al. [66] , 
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